Blood transcriptomics analysis offers insights into variant-specific immune response to SARS-CoV-2 (preprint)

Bulk RNA sequencing (RNA-seq) of blood is typically used for gene expression analysis in biomedical research but is still rarely used in clinical practice. In this study, we argue that RNA-seq should be considered a routine diagnostic tool, as it offers not only insights into aberrant gene expression and splicing but also delivers additional readouts on immune cell type composition as well as B-cell and T-cell receptor (BCR/TCR) repertoires. We demonstrate that RNA-seq offers vital insights into a patients immune status via integrative analysis of RNA-seq data from patients infected with various SARS-CoV-2 variants (in total 240 samples with up to 200 million reads sequencing depth). We compare the results of computational cell-type deconvolution methods (e.g., MCP-counter, xCell, EPIC, quanTIseq) to complete blood count data, the current gold standard in clinical practice. We observe varying levels of lymphocyte depletion and significant differences in neutrophil levels between SARS-CoV-2 variants. Additionally, we identify B and T cell receptor (BCR/TCR) sequences using the tools MiXCR and TRUST4 to show that - combined with sequence alignments and pBLAST - they could be used to classify a patients disease. Finally, we investigated the sequencing depth required for such analyses and concluded that 10 million reads per sample is sufficient. In conclusion, our study reveals that computational cell-type deconvolution and BCR/TCR methods using bulk RNA-seq analyses can supplement missing CBC data and offer insights into immune responses, disease severity, and pathogen-specific immunity, all achievable with a sequencing depth of 10 million reads per sample. Key PointsO_LIComputational deconvolution of transcriptomes can estimate immune cell abundances in SARS-CoV-2 patients, supplementing missing CBC data. C_LIO_LI10 million RNA sequencing reads per sample suffice for analyzing immune responses and disease severity, including BCR/TCR identification. C_LI

Gut microbial disruption in critically ill patients with COVID-19 associated pulmonary aspergillosis (preprint)

Objectives: COVID-19 disease can be exacerbated by Aspergillus superinfection (CAPA). The causes of CAPA are not yet fully understood. Recently, alterations in the gut microbiome have been associated with a complicating course and increasing severity of COVID-19 disease, most likely via immunological mechanisms. Aim of this study was to investigate a potential association between severe CAPA and alterations in the gut and bronchial microbiota. Methods: We performed 16S rRNA gene amplicon sequencing of stool and bronchial samples from a total of 16 COVID-19 patients with CAPA and 26 patients without CAPA. All patients were admitted to the intensive care unit. Results were carefully tested for potential influences on the microbiome during hospitalization. Results: We found that late in COVID-19 disease, CAPA patients exhibited a trend towards reduced gut microbial diversity. Furthermore, late stage CAPA disease showed an increased presence of Staphylococcus epidermidis in the gut. This is not found in late non-CAPA cases or early disease. The analysis of bronchial samples did not show significant results. Conclusions: This is the first study showing alterations in the gut microbiome accompany severe CAPA and possibly influence the hosts immunological response. In particular, an increase of Staphylococcus epidermidis in the intestine could be of importance.

Inference of differential gene regulatory networks from gene expression data using boosted differential trees (preprint)

Diseases can be caused by molecular perturbations that induce specific changes in regulatory interactions and their coordinated expression, also referred to as network rewiring. However, the detection of complex changes in regulatory connections remains a challenging task and would benefit from the development of novel non-parametric approaches. We developed a new ensemble method called BoostDiff (boosted differential regression trees) to infer a differential network discriminating between two conditions. BoostDiff builds an adaptively boosted (AdaBoost) ensemble of differential trees with respect to a target condition. To build the differential trees, we propose differential variance improvement as a novel splitting criterion. Variable importance measures derived from the resulting models are used to reflect changes in gene expression predictability and to build the output differential networks. BoostDiff outperforms existing differential network methods on simulated data evaluated in two different complexity settings. We then demonstrate the power of our approach when applied to real transcriptomics data in COVID-19 and Crohn’s disease. BoostDiff identifies context-specific networks that are enriched with genes of known disease-relevant pathways and complements standard differential expression analyses. BoostDiff is available at https://github.com/gihannagalindez/boostdiff_inference . Author Summary Gene regulatory networks, which comprise the collection of regulatory relationships between transcription factors and their target genes, are important for controlling various molecular processes. Diseases can induce perturbations in normal gene co-expression patterns in these networks. Detecting differentially co-expressed or rewired edges between disease and healthy biological states can be thus useful for investigating the link between specific disease-associated molecular alterations and phenotype. We developed BoostDiff (boosted differential trees), an ensemble method to derive differential networks between two biological contexts. Our approach applies a boosting scheme using differential trees as base learner. A differential tree is a new tree structure that is built from two expression datasets using a splitting criterion called the differential variance improvement. The resulting BoostDiff model learns the most differentially predictive features which are then used to build the directed differential networks. BoostDiff outperforms other differential network methods on simulated data and outputs more biologically meaningful results when evaluated on real transcriptomics datasets. BoostDiff can be applied to gene expression data to reveal new disease mechanisms or identify potential therapeutic targets.

Systematic analysis of alternative splicing in time course data using Spycone (preprint)

During disease progression or organism development, alternative splicing (AS) may lead to isoform switches (IS) that demonstrate similar temporal patterns and reflect the AS co-regulation of such genes. Tools for dynamic process analysis usually neglect AS. Here we propose Spycone ( https://github.com/yollct/spycone ), a splicing-aware framework for time course data analysis. Spycone exploits a novel IS detection algorithm and offers downstream analysis such as network and gene set enrichment. We demonstrate the performance of Spycone using simulated and real-world data of SARS-CoV-2 infection.

A Systematic Comparison of Differential Analysis Methods for CyTOF Data (preprint)

Cytometry techniques are widely used to discover cellular characteristics at single-cell resolution. Many data analysis methods for cytometry data focus solely on identifying subpopulations via clustering and testing for differential cell abundance. For differential expression analysis of markers between conditions, only few tools exist. These tools either reduce the data distribution to medians, discarding valuable information, or have underlying assumptions that may not hold for all expression patterns. Here, we systematically evaluated existing and novel approaches for differential expression analysis on real and simulated CyTOF data. We found that methods using median marker expressions compute fast and reliable results when the data is not strongly zero-inflated. Methods using all data detect changes in strongly zero-inflated markers, but partially suffer from overprediction or cannot handle big datasets. We present a new method, CyEMD, based on calculating the Earth Movers Distance between expression distributions that can handle strong zero-inflation without being too sensitive. Additionally, we developed CYANUS, a user-friendly R Shiny App allowing the user to analyze cytometry data with state-of-the-art tools, including well-performing methods from our comparison. A public web interface is available at https://exbio.wzw.tum.de/cyanus/.

Platelet expression and reactivity after BNT162b2 vaccine administration (preprint)

SARS-CoV-2 infection induces a coagulopathy characterized by platelet activation and a hypercoagulable state with an increased incidence of cardiovascular events. The viral spike protein S has been reported to enhance thrombosis formation, stimulate platelets to release pro-coagulant factors and promote the formation of platelet-leukocyte aggregates even in absence of the virus. Although SARS-CoV-2 vaccines induce spike protein overexpression to trigger SARS-CoV-2-specific immune protection, thrombocyte activity has not been investigated in this context. Here, we provide the first phenotypic platelet characterization of healthy human subjects undergoing BNT162b2 vaccination. Using mass cytometry, we analyzed the expression of constitutive transmembrane receptors, adhesion proteins and platelet activation markers in 12 healthy donors before and at five different timepoints within four weeks after the first BNT162b2 administration. We measured platelet reactivity by stimulating thrombocyte activation with thrombin receptor-activating peptide (TRAP). Activation marker expression (P-Selectin, LAMP-3, LAMP-1, CD40L and PAC-1) did not change after vaccination. All investigated constitutive transmembrane proteins showed similar expressions over time. Platelet reactivity was not altered after BNT162b2 administration. Activation marker expression was significantly lower compared to an independent cohort of mild symptomatic COVID-19 patients analyzed with the same platform. This study reveals that BNT162b2 administration does not alter platelet protein expression and reactivity.

Computational Strategies to Combat COVID-19: Useful Tools to Accelerate SARS-CoV-2 and Coronavirus Research (preprint)

SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) is a novel virus of the family Coronaviridae. The virus causes the infectious disease COVID-19. The biology of coronaviruses has been studied for many years. However, bioinformatics tools designed explicitly for SARS-CoV-2 have only recently been developed as a rapid reaction to the need for fast detection, understanding, and treatment of COVID-19. To control the ongoing COVID-19 pandemic, it is of utmost importance to get insight into the evolution and pathogenesis of the virus. In this review, we cover bioinformatics workflows and tools for the routine detection of SARS-CoV-2 infection, the reliable analysis of sequencing data, the tracking of the COVID-19 pandemic and evaluation of containment measures, the study of coronavirus evolution, the discovery of potential drug targets and development of therapeutic strategies. For each tool, we briefly describe its use case and how it advances research specifically for SARS-CoV-2. All tools are freely available online, either through web applications or public code repositories.

Exploring the SARS-CoV-2 virus-host-drug interactome for drug repurposing (preprint)

Coronavirus Disease-2019 (COVID-19) is an infectious disease caused by the SARS-CoV-2 virus. It was first identified in Wuhan, China, and has since spread causing a global pandemic. Various studies have been performed to understand the molecular mechanisms of viral infection for predicting drug repurposing candidates. However, such information is spread across many publications and it is very time-consuming to access, integrate, explore, and exploit. We developed CoVex, the first interactive online platform for SARS-CoV-2 and SARS-CoV-1 host interactome exploration and drug (target) identification. CoVex integrates 1) experimentally validated virus-human protein interactions, 2) human protein-protein interactions and 3) drug-target interactions. The web interface allows user-friendly visual exploration of the virus-host interactome and implements systems medicine algorithms for network-based prediction of drugs. Thus, CoVex is an important resource, not only to understand the molecular mechanisms involved in SARS-CoV-2 and SARS-CoV-1 pathogenicity, but also in clinical research for the identification and prioritization of candidate therapeutics. We apply CoVex to investigate recent hypotheses on a systems biology level and to systematically explore the molecular mechanisms driving the virus life cycle. Furthermore, we extract and discuss drug repurposing candidates involved in these mechanisms. CoVex renders COVID-19 drug research systems-medicine-ready by giving the scientific community direct access to network medicine algorithms integrating virus-host-drug interactions. It is available at https://exbio.wzw.tum.de/covex/.